Two-Phase Region of the Vortex-Sohd Melting Transition: 3D XY Theory 
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In clean enough samples of the high- Tc oxide materials, the phase transition into the super- 
conducting state occurs along a first order line in the H-T plane. This means that a two-phase 
region occurs in the B-T plane, in which the liquid and solid vortex phases coexist. We discuss the 
thermodynamics of this two-phase region, developing formulae relating experimental quantities of 
interest. We then apply the 3D XY scaling theory to the problem, obtaining detailed predictions for 
the boundaries of the coexistence region. By using published data, we are able to predict the width 
of the two-phase region, and determine the physical parameters involved in the 3D XY description. 
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The phase transition into the superconducting state in 
the high- Tc oxide materials is a topic of great current 
interest. It is well known that the mean-field descrip- 
tion of the superconducting transition must be modified 
to incorporate strong fluctuation effects. In the pres- 
ence of magnetic fields, the resulting transition has been 
described as the melting of a vortex-solid into a vortex- 
liquid, and is predicted to be first order in very clean 
samples p I . This prediction is borne out both in exper- 
iments [^5j and simulations [^. Important thermody- 
namic support for the first order scenario has been pro- 
vided by Schilling et ai, who verify that the Clausius- 
Clapeyron equation is satisfied 

In contrast, the zero field transition {H — Q,T = Tc) 
is continuous, and belongs to the universality class of the 
3D XY model The discontinuities of the entropy den- 
sity (As) and magnetization (AM) therefore both go to 
zero approaching this critical point. The observed range 
of the critical fluctuations extends over large portions of 
the phase diagram in both clean |^ and disordered ||] 
samples. The 3D XY theory provides a rigorous frame- 
work for analyzing the first order transition, and thus 
places strict constraints on its interpretation. Experi- 
mental data supporting this description have been ob- 
tained by Liang et al. 

As illustrated in Fig. 1, the first order transition occurs 
along a line Hm{T) in the H-T plane. However, in the 
B-T plane, the first order line becomes an area of two- 
phase coexistence. This area is analogous to the region of 
liquid-solid coexistence in the T-p (temperature-density) 
plane of a substance such as water. When B (the av- 
erage magnetic field) and T are such that the supercon- 
ductor is in the two-phase region, the system separates 
into two phases, each with its own value of B. It is a 
point of general and fundamental interest to understand 
the distinction between the H-T and B-T planes; as we 
shall explain, some nontrivial signatures of the first order 
transition can be inferred by keeping this distinction in 
mind. 

The purpose of this paper is to use 3D XY scaling 



to discuss the two-phase region. We start by giving a 
general thermodynamic analysis, emphasizing the Hc^-T 
plane, where H^^ is the externally applied magnetic field 
]lOt . This allows us to obtain relations between the dif- 
ferent experimental quantities of interest, thus providing 
important consistency checks for experimental investiga- 
tions. We then apply the 3D XY scaling theory to make 
more specific predictions concerning the coexistence re- 
gion. Background effects are carefully included, in terms 
of the background inverse permeability f2. We evalu- 
ate the important parameters of the scaling theory using 
published data; this allows us to make definite predic- 
tions about the size of the coexistence region for several 
systems. 

General Theory. We consider the coexistence of two 
phases: the vortex-liquid (L) and the vortex-solid (S). 
(These names are used only for convenience; the particu- 
lar nature of the two phases plays no role in the analysis.) 
We assume that each phase is characterized by a distinct 
free energy density, denoted as fL{B,T) and fs{B,T). 
The conjugate fields Ha{B,T) are defined in the usual 
way: 



Hc,{B,T)=47Tdfc./dB, 



(1) 



where a = L or S. The curves bounding the two- 
phase region in the B-T plane are denoted by B]^{T) 
and Bs{T), as in Fig. 1. The fundamental equations 
describing coexistence along the line H,n{T) are then 

Hl{Bl{T),T) = Hs{Bs{T),T) = H„,{T), (2) 
fs{Bs{T),T)~ fL{BL{T),T) 

= Hm{T)[Bs{T)-BL{T)]/A^. (3) 

The first equation requires that H be the same in the two 
phases at a given temperature, while the second equa- 
tion requires that the magnetic Gibbs free energies be 
the same at that temperature. Eqs. (^ and (|^) can be 
solved simultaneously to determine the two unknowns, 
Bs{T) and Bl{T). 

For most experimental situations, it is the external 
field ifex which is directly under control. For an ellip- 
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soidal sample with demagnetizing factor n, Hex is related 
to B and H by 

H^x = nB + {l-n)H = B + 4:TT{n- 1)M, (4) 

when the field is applied along a symmetry axis of the 
ellipsoid Q. 









\Hm(T) 






\bl(T) 


H 






Bs(T)\^ ^ 




T 


Tc 



Tc 

T 



FIG. 1. First order melting phase diagram in the B-T 
plane, showing the two-phase region (shaded). The coexis- 
tence boundaries are given as Bs{T) and Bl{T). Inset: in 
the H-T plane, coexistence occurs along the line Hm{T). 

We will now consider the two most common experi- 
mental procedures for collecting data: (i) varying Hex 
while keeping T fixed, (ii) varying T while keeping H^x 
fixed. It will be seen that the distinction between these 
two cases is nontrivial. 

The first results are evident from Eqs. ^ and (^) [p^ : 

AiJexIr = nAB|T = 47rnAAf |t. (5) 

Here, the operator A is defined as the value measured on 
the liquid side minus the value measured on the solid side 
of the transition. The notation of Eq. (||) also signifies 
that the difference is measured at constant T in this case. 
Note that the discontinuity AB\t occurs because fields 
in the range Bs{T) < B < Bl(T) are not stable. How- 
ever, the discontinuity AHg^ reflects only the boundaries 
of the two-phase region; there are no unattainable values 

of iJcx- 

We now discuss the constant-iJox path across the co- 
existence region. Eq. (|^) shows that a discontinuity 
AB\h^^ causes corresponding discontinuities in AH\h^^ 
and AM\h^^- According to Eq. (||), if iJ changes val- 
ues in the coexistence region, then it must follow the line 
Hm{T). Since the discontinuities AB and AH are much 
smaller than B and H respectively (except when T ~ Tc), 
we see that AH\h„^/ AT\h,^ dH„JdT. The following 
relations are thus obtained at constant Hex- 

AB|h„ = (1 - l/n)AH\H„^ = 47r(l - n)AM\H,^ 

~ (I - l/n)idH^/dT)AT\H^^. (6) 



In the final equation, H may be replaced by iJex using 
the approximate relation [1 -|- (n — 1)(1 — il,)]H ~ fi-ffcx, 
where the inverse background permeability ft is defined 
below, and satisfies O ~ I. Corrections to this approxi- 
mation are of order M/Hcx, and are therefore extremely 
small for typical YBa2Cu307_5 (YBCO) samples, but 
become more noticeable in Bi2Sr2CaCu208 (BSCCO). 

Equations (||) and (H) are related through the following 
geometrical statement: 

AH,x\t/ATU^ 2i -dH^xrrJdT. (7) 

This equation is important in its own right, since it is 
not affected by the sample geometry. 

Using the above results, we may obtain various rela- 
tions of interest. For example, we find [1 + (n — I)(l — 
n)]AM\H,^ =i flAMlr, with corrections of order M /Hex- 
When n — 1, this equation reduces to AM\h^^ — AM\t, 
which is independent of sample geometry. Similarly for 
n = 1, we find A{B — Hcx)\t — AB\h^^, which is also 
geometry- independent. Finally, we state the Clausius- 
Clapeyron equation appropriate for this transition: 

47rAs|T/AB|T = -dH„JdT. (8) 

Equations (|1)-(|1) form a set of thermodynamic relations 
between quantities of interest for the first order transi- 
tion. 

3D XY Scaling Theory. The main features of 3D XY 
scaling can be derived from a single ansatz for the free 
energy density [p^|-[T^. The new assumption for the co- 
existence region is that the liquid and solid free energies 
densities independently satisfy this scaling. This means 
that we may write the free energies in the following form, 
using t=\T~Te\/Te @: 

fL{B,T) - fkLt^''MBt-^''/HkL) + /o(T) + nBySn 
fs{B,T) = hst"'^s{Bt-^''/Hks) + /o(T) + rfS^/S^. 

Note the following points: 

(1) The exponent v ~ 0.67 is the correlation length 
exponent for the 3D XY model. The central assumption 
of the scaling ansatz is that the magnetic field B scales 
as an inverse length squared jl^. B and T then enter 
the scaling function only in the combination Bt^^^ . 

(2) The parameters Jul-, fks, HkL, and Hks are mate- 
rial dependent and thus nonuniversal. fk's have dimen- 
sions of free energy density, while H^s have dimensions 
of magnetic field. Below, we show that fkL = fks- 

(3) The term /o(T) + ^B'^/Stt approximates the 
smooth background contribution to the free energy den- 
sity near the zero field transition. The quantities foiT) 
and fl{T) are nonuniversal and contain no singularities 
near T ~Tc. The inverse background permeability fi,(T) 
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is only weakly temperature dependent in many cases, and 
satisfies ~ 1. We assume that the background terms 
have the same form in both the solid and liquid phases. 

(4) The two scaling functions (ps and 0^ are distinct 
and universal. However, the two phases should become 
indistinguishable when B — 0. Thus, we expect the fol- 
lowing equality to hold: 

fkLMO) = fks<PsiO). (9) 

(5) The conjugate fields are obtained from Eq. (|^): 
H^{B,T) ^nB + {ATrfkc./Hka.)f(bUBt-^-'/Hka) 

where a = L ot S. We now use Eqs. (1) and (|) to deter- 
mine the coexistence boundary lines Bs{T) and Bl{T). 
This is accomplished by noting that all terms in the 
free energy involving 0^ explicitly (the superconductiv- 
ity terms) are small compared to the remaining back- 
ground terms. (This assumption breaks down very near 
the zero field transition, where a more careful treatment 
is required [^.) The following explicit solutions are ob- 
tained: 

BLiT) = BH"^" + be, 
Bs{T) = B*t^'' -bt". 

According to our assumptions, B* b here, while the 
apparent negative value of Bs (T) (when T ~ Tc) is un- 
physical. The field-like parameters B* and b are deter- 
mined from the following equations: 

fks(t>s{B*/Hks) = fkLMB*/HkL), 
b - i27T/n)[{fks/Hks)^'s{B*/Hks) - 

UkLlHkL)4>'L{B*lHkL)]. 

The width of the coexistence region in the B-T plane 
can now be computed along the constant-T path: 

^B\T = BL{T)-Bs{T)=2bt\ (10) 

Similarly, along a constant-_B path, the temperature 
width of the coexistence region is given to leading order 

by 

AT|b = {Tcb/iyB*){B/B*)^^/^''^-^/\ (11) 

Using Eqs. (|1)-(||), these discontinuities can be related to 
other quantities of interest. For example, the magneti- 
zation discontinuity AM\t is found to scale in the same 
way as the magnetization in the 3D XY theory, as first 
pointed out by Liang et al. 

In the H-T plane, the phase transition is given by 

HmiT) = nBH^-' + [m+{AnhL/HkL)4>'L{B*lHkLW, 

where the term is again much smaller than the t^^ 
term. The coexistence boundaries in the H^^-T plane 



can then easily be computed using Eq. (^). The en- 
tropy discontinuity can be calculated to leading order 
from Eq. (^). For simplicity, we give the result in terms 
of B: 

As|t = (6B*n://7rT,)(B/B*)(3''-i)/2^. (12) 

Finally, we show that the material parameters fks and 
fkL must be equal. Eq. (||) can first be rewritten as 
fkhl fks = 0fcs(O)/0fcL(O). Since the scaling functions 
(pks and 4>kL are universal, the ratio 4>ks{^) / 't'khi^) must 
also be universal. fkL and fks are therefore related by 
a universal proportionality constant, which can be ab- 
sorbed into either or ipkL- Therefore fkL = fkS- 

Comparison with Experiments. We now use published 
data to obtain preliminary estimates for the 3D XY field 
parameters B* and b, and thus the width of the coexis- 
tence region, for several systems. Additionally, we check 
the thermodynamic relations (|l)-(0). In order to apply 
our analysis to real physical situations, it is necessary to 
assume a particular sample geometry. In accordance with 
the discussion above, we approximate the different sam- 
ples as ellipsoids, inscribed within the given dimensions. 
However, we note that several recent papers have consid- 
ered the differences between ellipsoidal geometries and 
more realistic plate-like geometries, which are more diffi- 
cult to treat exactly jlj,|l^ . Although geometrical effects 
of this type may partially account for the inconsistencies 
discussed below, we believe that some of the discrepancies 
may be too large to be attributable to geometry alone. 
We also make the approximation that fl = 1. 

The first experiments we shall consider, in very pure, 
optimally-doped YBCO single crystals, typically involve 
bulk magnetization or specific heat data obtained at con- 
stant T or iJex- In Ref. Liang et al. observe a dis- 
continuity in M which is well described by the 3D XY 
theory throughout the studied field range, 5-40 kG. (For 
simplicity, the different fields are all specified here in 
units of Gauss.) Using the results above, and assum- 
ing n ~ 0.93 as appropriate for this sample, we find that 
B* ~ 1.0 X 10^ G and 6 ~ 0.66 G. To take a particular 
example then, the width of the coexistence region should 
be equal to AT\h,^ =i 4 x 10"^ K when iJex = 40 kG. 
We have also obtained estimates for B* and b from two 
other published magnetization studies of YBCO, with re- 
sults nearly identical to those given above QJ^. These 
values may therefore be characteristic of optimally doped 
YBCO [11 . 

According to the estimates given above, the coexis- 
tence region may be difficult to resolve along the tem- 
perature axis, due to its narrow width. However in the 
experiments of Schilling et al. Q, estimates for AiJoxIr 
and AM It can be obtained simultaneously, thus allow- 
ing Eq. (||) to be tested. The following results are ob- 
served when T = 85 K: 4t:AM\t ~ 0.2 G and AiJoxIr ^ 
1500 G. (Similar results are obtained in Ref. Q. Note 
that the quantity AT studied in Ref. is unrelated 



3 



to the quantity AT\h^^ discussed above.) According to 
Eq. (||), the observed value of AH^xIt is too large by a 
factor of 7500. On the other hand, Schilling et al. have 
shown that the Clausius-Clapeyron equation is obeyed 
in this system. It seems likely therefore that the for- 
mer discrepancy does not signal the failure of the ther- 
modynamic relations; instead, it may only reflect sam- 
ple dependent effects, which may be generically identi- 
fied as "broadening." Further experiments are required 
to resolve this issue. In any event, the predicted width 
of the coexistence region in YBCO appears to be much 
narrower than the typical experimental resolution. Ob- 
served jumps in the different measured quantities should 
therefore be treated as true discontinuities. 

Bulk measurements of magnetization discontinuities 
have also been performed in pure BSCCO single crys- 
tals. (For example, Ref. fl^.) However, as our sec- 
ond example we consider the experiments of Zeldov et 
al., which involve small Hall sensors These devices 
are capable of measuring the local field B, thus avoiding 
many problems due to sample inhomogeneity. The fol- 
lowing approximate information can be obtained from 
Ref. |] for T ~ 80 K: AB|h„ ^ 0.4 G, AT\h,^ < 
3 mK, A{B - Hex) It ^ 0.4 G, AiJexIr < 0.4 G, and 
dHcxm/dT ~ —6 G/K. The appropriate demagnetiz- 
ing factor for this sample is n ~ 0.71. Using these 
estimates, the geometry-independent (51 ~ 1) relation 
A{B — i?cx)|T — ^B\h^^ is well satisfied. The relation 
A{B — Hox)\t — — (1 — l/"-)AiJcx|T is not quite satisfied; 
this small discrepancy may be attributable to geomet- 
ric effects. However, relation (||) between AB\h^^ and 
AT\h^^ shows that the observed discontinuity AT\h^^ is 
at least 50 times sharper than predicted. Similarly, if the 
estimates given above for AT\h^^ and ATJcxIt are taken 
as equalities, the geometry-independent result of Eq. (|^) 
is violated by a factor of more than 20. To improve this 
situation, we would need to assume that AHcx\t was 
20 times smaller than estimated. Unfortunately, such a 
modification leads to discrepancies in other thermody- 
namic relations. 

The discrepancies observed in the local B measure- 
ments of BSCCO are of a different nature than those in 
YBCO. In the latter case, the observed width of the two- 
phase region was much too large, a fact which can be at- 
tributed to broadening. In the BSCCO case, the observed 
two-phase region was much too narrow, a result which 
cannot be explained by any obvious thermodynamic ar- 
guments. However, it is instructive to continue with the 
present analysis to obtain the 3D XY parameters appro- 
priate for BSCCO. According to Ref. Q, a reasonable 
estimate for B* is given by B* ~ 1000 G. This result is 
not affected by the inconsistencies described above, and is 
corroborated by other experiments. Assuming the value 
given above for A{B — -ffox)|T to be the most reliable es- 
timate of the discontinuities in BSCCO, we then obtain 
b ~ 2.8 G. From Ref. iQ (using low transverse fields. 



Fig. 10) we obtain a slightly different result of 6 = 0.8 G. 
Conclusions. We conclude by discussing two points. 

(1) The results given above for first order melting near 
the 3D XY critical point differ in several ways from ear- 
lier analyses aimed at fluid- like critical points p^ ]: (i) in 
the 3D XY scaling ansatz, it is the density-like variable 
(_B), not the pressure-like variable (H), which appears 
naturally in the scaling functions; (ii) in contrast with 
the fluid system, neither of the scaling variables B or H 
are directly related to the superconducting order param- 
eter; (iii) although a relation like Hks — HkL occurs in 
fluid- like scaling theories, the same cannot be proven for 
the vortex-solid melting transition in the absence of a 
more microscopic theory. 

(2) Our analysis should also apply to computer simu- 
lations aimed at understanding the first order transition 
Q . Many of these simulations are performed by varying 
the temperature at fixed B. In this two-phase re- 
gion should be observed, with width AT\b. If the value 
of B is small enough for 3D XY scaling to be valid, this 
width should be described by Eq. (pi]), and the entropy 
jmnp by Eq. (H). 

In summary, we have presented both general and spe- 
cific (3D XY) theories of the two-phase region for the 
vortex-solid melting transition. We find that in typical 
YBCO samples, the solid-liquid coexistence region may 
be too small to be observed experimentally. However, 
based upon available estimates for the discontinuity of 
B, the region should be easily detected in BSCCO using 
precision measurements. The observations of two-phase 
coexistence and thermodynamic self-consistency are both 
crucial for confirming the first order nature of the melt- 
ing transition. The recent experiment of Schilling et al. 
provides a check of the self-consistency relations; in this 
case the Clausius-Clapeyron equation is well satisfied, 
although open questions remain, regarding other rela- 
tions. It is hoped that future experiments may address 
this issue. Additionally, the scaling of the coexistence re- 
gion provides information regarding the most fundamen- 
tal fluctuations associated with superconductivity, which 
are governed by the 3D XY critical point. 
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